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Abstract 

The spectrum of extra-galactic cosmic rays (CRs) is expected to follow the Greisen- 
Zatsepin-Kuzmin (GZK) cutoff at about 5 x 10 10 GeV which results from energy losses of 
charged nuclei in the cosmic microwave background. So far the confrontation of this feature 
with CR data is inconclusive. In the absence of close-by sources a power-law continuation 
of the spectrum might signal the contribution of new physics. We have investigated the 
statistical significance of a model where exotic interactions of cosmogenic neutrinos are 
the origin of super-GZK events. A strong neutrino-nucleon interaction is favored by CR 
data, even if we account for a systematic shift in energy calibration. 



1. Introduction 

The appearance of extremely high energetic (EHE) cosmic rays (CRs) is a mys- 
tery. Simple geometric arguments show that the energy of nuclei originating in cos- 
mic accelerators should be limited corresponding to the size and magnetic field of 
the accelerating environment (Hillas criterion) [1]. Additional constraints arise from 
energy losses in the source that has to be balanced by the acceleration rate [2] . In 
particular, the energy loss length of synchrotron radiation decreases with the third 
power of the magnetic field (at constant gyro-radius) and limits the efficiency of 
small candidate acceleration sites. The very few sources that seem to be capable of 
accelerating protons up to 10 12 GeV include radio galaxy lobes and relativistically 
moving sources like jets of active galactic nuclei or gamma ray bursts. 

Not only the acceleration mechanism of these particles seems to be problem- 
atic, but also their propagation in the interstellar medium. Resonant pion photo- 
production in the cosmic microwave background (CMB) limits the range of EHE 
particles to a few 10 Mpc. It was first emphasized by Greisen [3], Zatsepin, and 
Kuzmin [4] (GZK) that this should be signaled by a cutoff in the CR spectrum at 
about 5 x 10 10 GeV if extra-galactic sources dominate at these energies. So far, the 
confrontation of this effective cutoff with CR data is inconclusive. In particular, the 
observations of two large exposure experiments AGASA [5] and HiRes [6] seem to 
imply conflicting results. 

The "excess" of super-GZK events reported by the AGASA collaboration has 
led to speculations about a different origin of EHE CRs. Berezinsky and Zatsepin 
proposed that cosmogenic neutrinos [9] produced in the decay of the GZK pions 
could explain these events assuming a rapid rise of the neutrino-nucleon interaction. 
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Fig. 1. Left panel: The cosmic ray spectra from AGASA, Fly's Eye Stereo, Hires-I/II (with energy 
shift given in Ref. [14]), and PAO. Right panel: The flux of extra-galactic protons and cosmogenic 
neutrinos from CMB interactions corresponding to the 99% CL of the fit. 

We have followed this idea in Ref. [11] and investigated the statistical significance 
of scenarios with strongly interacting neutrinos combining CR data from AGASA 
and HiRes and neutrino limits from horizontal events at AGASA and contained 
events at RICE. 

The AGASA excess of super-GZK events could also be the result of a relative 
systematic error in energy calibration of ±30% [12]. This interpretation is fueled 
by a recent re-analysis of the AGASA data [13]. The preliminary results indicate 
that the energy of highly inclined showers has previously been over-estimated. If 
verified, this would result in a re-calibration of the spectrum according to a relative 
energy shift of -10% to -15%. 

Remarkably, this energy shift agrees with predictions from a matching procedure 
of CR spectra assuming an early onset in extra-galactic proton dominance at about 
10 9 GeV [14]. The calibration by the "dip" in these proton spectra resulting from 
electron-positron pair production on the CMB predicts a —10% and 20% energy 
shift of the AGASA and HiRes data, respectively. We have incorporated this re- 
calibration into our analysis and show also the effect of the preliminary data from 
the Pierre Auger Observatory (PAO) [15] on our results. 



2. Cosmic Sources of Protons and Neutrinos 

We assume that extra-galactic protons dominate the CR spectrum above 5 x 
10 s GeV. A convenient parameterization of a spatially homogeneous and isotropic 
source luminosity is given by 

£ G r(z,E) oc {1 + z) n e"^, z min < z < z max , (1) 

which accounts for an evolution with redshift z and an exclusion of nearby ( z m i n ) 
and early (z max ) sources. In our analysis we have keep these parameters fixed at 
Zmin = 0.012, corresponding to r m ; n w 50Mpc, and z max = 2.0. The power-law 
injection is effectively limited to energies below E max , which we also fix at 10 12 GeV. 

The flux of protons originating at distant sources is subject to energy redshift and 
collisions with the interstellar photon background during propagation. The domi- 
nating interactions of UHE protons are e + e~ pair production and meson photopro- 
duction in the CMB (see e.g. [16]). These effects are taken into account by means of 



2 



propagation functions, which have been provided by the authors of Ref. [21]. The 
pions, which are resonantly produced at the GZK cut-off, decay into electron and 
muon neutrinos via the reaction chain n + — > n + v^ — > v^v^v e e + and the conjugate 
process. This constitutes a flux of UHE cosmogenic neutrinos which can be taken as 
a "guaranteed source" assuming extra-galactic protons as the highest energy CRs. 

Depending on the environment of the proton accelerators neutrinos might also 
emerge in the decay of pions produced by photoproduction processes during accel- 
eration [23]. The flux of extra-galactic neutrinos have been estimated in Ref. [24] 
assuming a low cross-over to extra-galactic protons at about 5 x 10 8 GeV. In con- 
trast to a transition at the "ankle" this flux of neutrinos from the source dominate 
the total neutrino flux below 10 9 GeV and exhausts current experimental neutrino 
bounds at energies of about 10 7 GeV. The total neutrino flux beyond the GZK 
cutoff is however dominated by cosmogenic neutrinos and we will only use this 
contribution for our statistical evaluation. 



3. Exotic Neutrino Interactions 

The flux of cosmogenic neutrinos associated with a low cross-over to extra-galactic 
protons is comparable to the flux of protons at the GZK cut-off as can be seen from 
Fig. 1. A contribution of these neutrinos as super-GZK events requires a strong 
deviation from the feeble neutrino-nucleon interactions predicted by the Standard 
Model. In order to avoid large contributions in neutrino experiments the transition 
to a anomalously large cross section has to be very rapid. 

The realization of such a behavior has been proposed in scenarios beyond the 
(perturbative) Standard Model (SM), e.g. arising through compositeness, through 
electroweak sphalerons, through string excitations in theories with a low string and 
unification scale, through Kaluza-Klein modes from compactified extra dimensions, 
or through black hole and p-brane production, respectively. For details of these 
models we refer to some recent reviews, Refs. [25]. Skepticism about these scenarios 
has been raised in Ref. [27] based on limitations from s-wave unitarity or the natu- 
ralness of the couplings. However, these considerations do not necessarily apply for 
non-perturbative aspects like instantons. 

In the following we will use a flexible parameterization of a strong neutrino- 
nucleon inelastic cross section (ct"™) focusing on three characteristic parameters: 
(i) the energy scale E t h of the new underlying physics, (ii) the amplification A 
compared to the SM predictions, and (iii) the width B of the transition between 
weak and strong interaction. A mathematical convenient parameterization is given 
by 



log 



' ' A*™) ~ 2 



E 

1 + tanh ( log e 



(2) 



In general, experiments distinguish different CR primaries by their characteristic 
shower development in matter. For the sake of simplicity we will assume that the 
characteristics of the showers induced by strongly interacting neutrinos are indis- 
tinguishable from those induced by protons. In particular, we assume for both 
primaries (i) a complete conversion of the incident energy into the shower, and (ii) 
equal detection efficiencies at the highest energies. 
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Fig. 2. The range of the cross section within the 95%, 90% and 68% CL (outer to inner bands) ob- 
tained by the re-calibrated AGASA, Fly's Eye Stereo, and HiRes-I/II data, with (right panel) and 
without (left panel) the preliminary Auger data. Also shown are theoretical predictions of the neu- 
trino-nucleon cross section enhanced by electroweak sphalerons, p-branes, and string excitations 
(see Ref. [11] for details). 



4. Quantitative Analysis 

Following the procedure of Ref. [11] we have evaluated the CR spectra of AGASA, 
Fly's Eye Stereo, HiRes-I/II, and PAO (preliminary data) assuming extra-galactic 
proton dominance above and 5 x 10 8 GeV and a contribution of strongly interacting 
neutrinos as super-GZK events. In contrast to Ref. [11] we have re-calibrated the 
energy by a multiplicative factor 0.9 and 1.2 for the AGASA and HiRes data, 
respectively, following Ref. [14]. For consistency with neutrino experiments we have 
used the search results of horizontal events at AGASA and contained events at 
RICE [28]. Figure 2 shows the range of the cross section corresponding to the 68%, 
90%, and 95% confidence level (CL) from a goodness-of-fit test assuming cosmogenic 
neutrino fluxes from a low crossover scenario. The details of the statistical analysis 
and the approximations involved can be found in Ref. [11]. 

Compared to our previous results [11] the re-calibrated spectrum is in a much 
better statistical agreement with the combined flux of protons and neutrinos from 
exotic interactions. This is mainly due to the fact that the calibration method 
"optimizes" the CR data to extra-galactic proton spectra below the GZK cutoff. 
However, within this particular model the re-calibrated data still favors an addi- 
tional contribution from strongly interacting cosmogenic neutrinos. At the 95% CL 
the corresponding neutrino-nucleon cross section should exhibit a steep increase by 
an amplification factor of A > 10 3 . For A < 10 5 the transition should be very rapid 
(B < 10) at about 5 x 10 11 GeV. The allowed parameter space further shrinks if 
we also include the preliminary PAO data. 

As an illustration of the result, we have considered three models of a rapidly 
increasing neutrino-nucleon cross section based on electroweak sphalerons [29], p- 
branes [30] and string excitations [31]. The details can be found in Ref. [11]. A 
separate fit of the source luminosity (Eq. (1)) with these cross section using the re- 
calibrated data gives a statistical acceptance at the 81%, 65%, and 60% CL (without 
PAO), respectively. If we also include the PAO data in the fit, these values change 
to 98%, 90%, and 88%, respectively, which is in reasonable good agreement with 
the results of Fig. 2. 
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5. Conclusions 

We have shown that in the absence of close-by sources exotic interactions of 
cosmogenic neutrinos might extend the cosmic ray spectrum beyond the GZK cut- 
off. Such a behavior is predicted in various extensions of the Standard Model, 
but might also be due to non-perturbative aspects like electro-weak sphalerons. In 
extension to our analysis of Ref. [11] we have presented here the results of a fit to 
the AG AS A and HiRes data, shifted according to the "dip" -calibration from extra- 
galactic proton spectra [14]. With these modifications a contribution of neutrinos 
is still favored by the data. The inclusion of preliminary data of the Pierre Auger 
Observatory weakens the results of the goodness-of-fit test. 
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